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PARAMETRIC  INVESTIGATION  OF  RADOME  ANALYSIS  METHODS 

I.  Introduction  and  Summary 

The  overall  objective  of  this  research  is  to  develop  a  general 
theory  of  radome  analysis  and  to  determine  the  accuracies  of  various 
radome  analysis  methods  under  controlled  conditions  of  antenna  size 
and  placement,  wavelength,  and  radome  size  and  shape.  Experimental 
measurements  on  selected  antenna/radome  combinations  at  35  GHz  will  be 
used  as  true  data  in  the  assessment  of  accuracy. 

During  the  first  year  a  general  theory  of  analysis  was  developed 
based  on  the  Lorentz  reciprocity  theorem  and  on  the  Huygens-Fresnel 
principle.  Computed-aided  methods  of  analysis  based  on  these  principle 
were  formulated  and  programmed.  Three  monopulse  antennas  were  designed 
and  fabricated  for  use  in  carrying  out  the  measurements  at  35  GHz. 

Three  radomes  were  also  designed  and  fabricated  for  use  in  the  experi¬ 
mental  program.  A  mechanical  fixture  was  fabricated  which  allows  accu¬ 
rate  positioning  of  the  antenna  with  respect  to  the  radome  and  with 
respect  to  the  reference  system  used  for  pattern  measurements.  A 
detailed  discussion  of  the  progress  in  the  first  year  was  presented  in 
Reference  1. 

During  the  second  year,  the  initiation  efforts  of  the  first  year 
work  have  begun  to  come  to  fruition  through  concentrated  experimental 
and  analytical  efforts  to  gather  true  and  calculated  data  which  will 
serve  as  the  basis  for  realizing  the  overall  research  objective.  Ex¬ 
tensive  pattern  measurements  have  been  made  for  nine  of  the  fifteen 


antenna/radome  combinations  available.  These  experimental  efforts  have 


been  paralleled  by  computed-aided  analyses  of  the  same  combinations 
for  comparison  purposes.  In  the  process,  a  technique  for  synthesizing 
the  aperture  fields  of  the  four-horn  monopulse  antennas  from  measured 
principal  plane  patterns  has  been  developed.  The  progress  during  the 
second  year  is  described  in  Section  II  below. 

In  the  third  year,  measurements  will  be  carried  out  on  lens 
radomes  and  on  conventional  radomes  to  which  metal  tips  have  been  added, 
and  comparisons  to  theoretical  predictions  will  be  made.  Antenna/radome 
configurations  which  lend  themselves  to  exact  analysis  using  boundary¬ 
valued  approaches  and  modal  expansions  will  be  investigated. 

Four  publications  and  three  presentations  concerning  this  re¬ 
search  have  resulted.  They  are  described  in  Section  III  along  with 
papers  planned  for  submission  to  technical  journals. 

The  professional  personnel  associated  with  the  research  effort 
are  listed  in  Section  IV  along  with  the  iteractions  with  DOD  agencies 
and/or  their  contractors. 

Appendix  A  contains  copies  of  the  papers  presented  at  symposia. 
Appendix  B  contains  Chapter  1  of  a  planned  scientific  report  to  be 
submitted  for  review  in  March. 

II .  Progress  to  Date 

Efforts  during  the  second  year  have  been  made  in  both  analytical  and 
experimental  areas.  In  the  analytical  area,  the  radome  analysis  computer 
programs  have  been  finalized,  and  a  method  of  monopulse  antenna  snythesis 
has  been  implemented  to  greatly  facilitate  the  comparisons  of  computed  and 
neasured  effects  of  the  radomes  on  the  antenna  radiation  patterns.  In  the 
experimental  area,  radiation  pattern  and  boresight  error  measurements  have 
been  completed  on  the  various  combinations  of  antennas  and  radomes. 


The  radome  analysis  based  on  the  Huygens-Fresnel  principle  [2] 
and  Lorentz  reciprocity  [3]  has  been  fully  implemented  in  software  and 
successfully  tested  for  accuracy.  The  inner  radome  surface  is  used  as 
the  surface  of  integration.  The  normal  voltage  transmission  coeffi- 
ceints  are  used  to  transform  the  incident  plane  wave  from  the  outside 
surface  to  the  inside  surface  s'  where  the  surface  integration  is  per¬ 
formed  to  obtain  the  voltage  received  by  the  antenna  according  to 
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where  the  prime  indicates  integration  over  the  radome  surface. 

Testing  of  the  surface  integration  algorithms  has  been  done  by 

using  an  infinitely  thin,  "air"  radome  so  that  the  surface  integration 

is  carried  out  over  a  tangent  ogive  surface  in  free  space.  The  antenna 

is  located  at  the  base  of  the  radome  and  consists  of  a  square  aperture 

excited  with  a  TE,„  waveguide  mode  so  that  the  aperture  fields,  E  ,  H  , 
10  -ap  -ap 

are  approximately  Maxwellian.  A  plane  wave  is  then  assumed  incident 
on  the  antenna/radome  combination  from  the  direction  k.  The  received 
voltage  is  computed  according  to  Equation  (1)  and  according  to 
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where  S  represents  the  aperture  surface.  By  the  Lorentz  reciprocity 


theorem. 
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Table  1  shows  some  results  of  this  testing.  A  hemisphere  (F=.5) 
and  a  tangent  ogive  surface  with  fineness  ratio  F=1.0  were  used.  The 
angles  <fi  and  0  specify  the  direction  of  arrival  of  the  vertically  pola¬ 
rized  plane  wave  on  the  antenna  (z  normal  to  the  y-polarized  aperture) . 
Agreement  is  excellent. 

The  fields  E^,  used  in  Equation  (1)  were  computed  from  the 

specified  aperture  fields  using  the  Huygens-Fresnel  principle.  Both 

magnetic  (E  X  n)  and  electric  (n  X  H  )  currents  are  used.  Initially, 
— ap  -ap 

it  was  intended  that  the  plane  wave  spectrum  (PWS)  formulation  [4]  be 
used  to  compute  these  fields;  however,  its  implementation  has  not  yet 
been  successful  thus  far.  It  has  been  established  in  a  mathematically 
rigorous  manner  that  the  two  formulations  are  entirely  equivalent  and 
should  indeed  yield  the  same  results.  Work  is  continuing  in  this  area 
because  the  PWS  approach  may  offer  a  speed  advantage  in  the  computations. 

Another  important  analytical  endeavor,  and  one  that  has  required 
considerable  time,  is  the  development  and  implementation  of  an  aperture 
synthesis  method  using  measured  principal  plane  patterns  as  described 
in  Appendix  A.  Using  this  technique,  antenna  aperture  fields  are  gene¬ 
rated  from  the  measured  far-field  patterns  such  that  when  they  are  used 
in  the  radome  analysis,  the  measured  patterns  will  be  predicted  exactly 
in  the  absence  of  the  radome.  Then,  when  the  patterns  are  computed 
with  the  radome  in  place  and  compared  to  their  experimental  counterparts, 
valid  comparisons  can  be  made  between  the  accuracies  of  three  radome 
analysis  methods  under  investigation.  The  synthesis  technique  has  much 
broader  application  and  is  currently  being  investigated  for  two-dimen¬ 
sional  synthesis  using  measured  near-field  data  [5]  . 

Numerous  pattern  and  boresight  error  measurements  have  been  made 

on  the  eight  antenna/radome  combinations  shown  in  Table  2.  (The  large 
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Table  1.  Comparisons  of  Results  Obtained  for 
Radome  Surface  Integration  Testing 


Finess 

Ratio 

(deg. ) 

0 

(deg. ) 

V 

Amp. (dB) 

t 

R  Q 

Phase (  ) 

Amp(dB) 

V 

R 

Phase (  ) 

No.  of  Points/ 
Computation  Timet  sec.) 

.5 

0 

0 

0.00 

180. 

0.00 

180. 

162/11.69 

.5 

0 

15 

-4.94 

179. 

-5.11 

-180. 

162/12.01 

.  5 

90 

15 

-8.22 

-180. 

-8.32 

-180. 

162/12.33 

1.0 

0 

0 

0.00 

180. 

0.00 

180. 

229/17.53 

1.0 

0 

15 

-4.95 

180. 

-5.11 

-180. 

229/17.91 

1.0 

90 

15 

-8.25 

-179. 

-8.32 

-180. 

229/17.48 

Table  2.  Ratios  of  D  /D  for  Selected  Antennas  and  Radomes 

is  ap 


Antenna 


Radome 

Smali 

Medium 

Large 

Small  ( F=l) 

2.33* 

— 

-- 

Medium  (F=l) 

3.98* 

2.33* 

1.27 

Medium  (F=1.5) 

— 

2.33* 

— 

Medium  (F=2.0) 

— 

2.33* 

— 

Large  (F=l) 

7.28* 

4.27* 

2.33* 

*  Measurements  completed. 


antenna  would  not  fit  inside  the  medium  radome.)  For  each  antenna/ 
radome  combination,  eight  patterns  were  measured  for  each  of  the  three 
channels  of  the  four-horn  monopulse  antenna:  four  pattern  cuts  (4>=0, 

45°,  90°,  and  135°)  and  two  polarizations.  In  addition,  boresight  error 
measurements  were  made  for  each  combination  as  allowed  by  the  special 
positioner  used. 

All  the  patterns  were  made  on  a  pattern  range  in  the  usual  manner; 
however,  considerable  effort  was  spent  in  translating  the  data  from 
pattern  recording  paper  to  computer  files  for  easy  access  during  the  com¬ 
parison  phase.  This  translation  was  done  using  an  xy  digitizer  and  HP2100 
minicomputer  system  at  the  Engineering  Experiment  Station.  After  the 
data  were  digitized  and  stored  in  files  on  the  HP2100,  they  were  trans¬ 
ferred  to  the  central  campus  computer  (Cyber  70/74)  where  the  analyses 
are  being  done. 

Plotting  software  has  been  developed  on  the  Cyber  for  overlaying 
various  patterns  on  the  same  axes  for  comparison  purposes .  An  example 
is  shown  in  Figure  1  for  the  azimuth  difference  channel  of  the  medium 
antenna  inside  the  medium  radome:  the  solid  line  is  the  pattern  with 
the  radome;  the  dashed  is  the  pattern  of  the  antenna  without  the  radome 
(for  reference) . 

At  the  present  time,  the  aperture  fields  for  the  three  antennas 
are  being  synthesized  using  measured  pattern  data  such  as  that  shown  in 
Figure  1.  The  analyses  will  be  run  next  and  comparisons  made  between 
measured  and  computed  results.  From  these  comparisons,  the  accuracies 
and  ranges  of  validity  of  the  fast  receiving,  fast  transmitting,  and 
surface  integration  methods  of  analysis  will  be  determined. 

During  the  third  year  (which  is  already  in  progress) ,  three  topics 
are  of  interest:  radomes  with  metal  tips,  lens  radomes,  and  antenna/ 


radome  configurations  which  lend  themselves  to  exact  analysis.  The 


three  medium  sized  radomes  having  fineness  ratios  of  1,  1.5,  and  2  will 
be  fitted  with  metal  tips,  and  measurements  and  analytical  comparisons 
will  be  made.  The  analysis  software  will  be  modified  to  accommodate 
lens  radomes,  and  measurements  will  be  carried  out  on  two  lenses.  Finally 
the  boundary- valve  problem  of  a  perfectly  conducting  sphere  inside  a  con¬ 
centric  spherical  dielectric  shell  will  be  solved  to  serve  as  a  model 
for  generating  exact  analytical  data  for  general  use  in  accessing  the 
accuracy  of  any  computer-aided  radome  analysis  method. 

Ill .  Presentations  and  Publications 


G.  K.  Huddleston,  H.  L.  Bassett,  and  J.  M.  Newton,  "Parametric 
Investigation  of  Radome  Analysis  Methods,"  Presented  at  and 
published  in  Proceedings  of  1978  IEEE  AP-S  Symposium  and 
Proceedings  of  Fourteenth  Symposium  on  Electromagnetic  Windows. 

G.  K.  Huddleston,  "Radome  Analysis,"  invited  presentation  to 
Atlanta  Chapter  of  IEEE,  AP-S/MTT  Group. 

G.  K.  Huddleston,  H.  L.  Bassett,  and  J.  M.  Newton,  "Parametric 
Investigation  of  Radome  Analysis  Methods,"  Annual  Technical 
Report,  AFOSR*  November  1978. 

G.  K.  Huddleston,  "Equivalence  of  Receiving  and  Transmitting 
Formulations  of  Radome  Analysis,"  Paper  to  be  submitted  to 
IEEE  AP-S  Transactions  (see  Appendix  A) . 

G.  K.  Huddleston,  "A  General  Theory  of  Radome  Analysis,"  Paper 
in  preparation  for  submission  to  IEEE  AP-S  Transactions. 

G.  K.  Huddleston  and  I.  Rusodimos,  "Computed- Aided  Radome  Analysis 
Using  Geometrical  Optics  and  Lorentz  Reciprocity,"  AFOSR  Scien¬ 
tific  Report,  planned  for  submission  in  March  1980  (see  Appendix 
B)  . 


G.  K.  Huddleston  and  I.  Rusodimos,  "Computed-Aided  Radome  Analysis 
Using  Huygens-Fresnel  Principle  and  Lorentz  Reciprocity,"  AFOSR 
Scientific  Report,  planned  for  submission  in  April  1980. 

G.  K.  Huddleston,  "Aperture  Synthesis  From  Principal  Plane  Patterns 
of  Monopulse  Antenna  for  Radome  Analysis,"  submitted  to  1980 
IEEE/AP-S  Symposium,  June  1980  (see  Appendix  A) . 


G.  K.  Huddleston,  H.  L.  Bassett,  and  J.  M.  Newton,  "Parametric 
Investigation  of  Radome  Analysis  Methods:  Experimental  Results, 
AFOSR  Scientific  Report,  planned  for  submission  April  1980. 


G.  K.  Huddleston,  H.  L.  Bassett,  and  J.  M.  Newton,  "Parametric 
Investigation  of  Radome  Analysis  Methods:  Comparisons  of 
Analytical  and  Experimental  Results,"  AFOSR  Scientific  Report, 
planned  for  submission  June  1980;  also  planned  for  submission 
as  paper  to  IEEE  Transactions  on  Antennas  and  Propagation. 

IV.  Personnel  and  Interactions 

The  following  professional  personnel  are  actively  engaged  in 

this  research  program: 

G.  K.  Huddleston  Assistant  Professor  School  of  Electrical 

Engineering 

H.  L.  Bassett  Senior  Research  Engineering  Experi- 

Engineer  ment  Station 

J.  M.  Newton  Research  Engineer  Engineering  Experi¬ 

ment  Station 

Dr.  Huddleston  and  Mr.  Bassett  serve  as  co-principal  investigators. 
Dr.  Huddleston  directs  the  analytical  efforts  while  Mr.  Bassett  directs 
the  experimental  work.  Mr.  Jason  Rusodimos,  a  graduate  student  in  the 
School  of  Electrical  Engineering  worked  with  Dr.  Huddleston  on  the  com¬ 
puter  implementation  of  various  analysis  methods.  Mr.  William  Adams, 
a  co-operative  student,  has  worked  closely  with  Mr.  Newton  in  the  mea¬ 
surements  program. 

Interactions  have  taken  place  with  DOD  agencies  and/or  their 
contractors.  Software  has  been  delivered  to  Naval  Surface  Weapons  Center, 
Dahlgren,  Virginia,  (Richard  M.  Swink)  for  use  in  assessing  the  performance 
of  a  new  shipboard  radome  design.  Radome  analysis  computer  programs 
have  been  documented  and  installed  on  the  IBM  3033  system  at  Applied 
Physics  Laboratory,  Laurel,  Maryland,  (R.  C.  Mallalieu)  for  use  in  a 
significant  radome  technology  program  for  the  Department  of  the  Navy.  The 
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fast  receiving  formulation  software  has  been  delivered  to  Flight  Systems, 
Inc.,  Newport  Beach,  California,  (Bernard  Crowe)  for  general  radome 
analysis  work  in  the  ECM  area.  This  software  has  also  been  installed 
on  the  SEL  system  at  Georgia  Tech  Engineering  Experiment  Station  (H.  L. 
Bassett)  and  is  being  used  in  a  polarimetric  study  for  the  U.  S.  Army 
Missile  Command.  Discussions  are  currently  underway  with  China  Lake 
Test  Center  (Joseph  Mosko)  and  Texas  Instruments,  Inc.  (Dave  Burks) 
concerning  this  radome  work. 

Three  other  developments  are  noteworthy.  First,  the  School  of 
Electrical  Engineering  has  recently  become  a  member  of  the  Joint  Services 
Electronics  Program,  and  Dr.  Huddleston  will  be  carrying  out  an  investiga¬ 
tion  to  determine,  by  measurement,  the  electromagnetic  fields  on  the 
surface  of  a  radome.  Second,  a  specialized  research  equipment  grant 
has  been  received  from  the  National  Science  Foundation  to  establish  an 
Automated  Radiation  Laboratory  in  the  School  of  Electrical  Engineering 
for  use  in  carrying  out  antenna  and  radome  research.  Finally,  an 
advanced  graduate  course  on  radome  analysis  was  taught  last  Fall  by 
Dr.  Huddleston  (see  Appendix  C) . 

It  is  important  to  note  that  the  AFOSR  grant  support  has  been 
catalytic  and  instrumental  in  the  success  of  these  auxiliary  radome- 
related  activities,  and  that  support  is  hereby  gratefully  acknowledged. 
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Valid  comparisons  of  the  accuracies  of  radorae  analysis  [l]  using 
measured  radiation  patterns  as  true  data  require  that  each  analysis 
method  be  capable  of  predicting  exactly  the  antenna  patterns  measured 
in  the  absence  of  the  radome.  When  the  measured  amplitude  and  phase 
of  the  antenna  fields  are  known  over  the  complete  sphere,  an  aperture 
synthesis  method  such  as  that  described  by  Ludwig  [2]  can  be  used. 

But  when  the  measured  data  is  limited  to  amplitude-only,  principal 
plane  patterns,  a  different  approach  and  one  which  exploits  all 
available  information  about  the  antenna  to  be  synthesized,  must  be 
used. 


Such  a  method  of  aperture  synthesis  for  a  four-horn  monopulse 
antenna  using  principal  plane  amplitude  patterns  of  the  sum  and  two 
difference  channels  is  described.  The  radius  (a11. 74 A)  and  element 
half-spacing  (dx=dy=.95A)  of  each  vertically  (y)  polarized  conical 
horn  element  (10°  flare  angle)  are  also  used  in  the  procedure.  The 
accuracy  of  the  method  is  illustrated  below  for  a  theoretical  array. 
Synthesis  results  obtained  using  measured  data,  and  analysis  results 
with  the  radome  in  place,  will  be  presented  at  the  symposium. 

Theory 

The  synthesis  procedure  is  carried  out  independently  in  the  two 
principal  planes  to  determine  the  plane  wave  spectra  Axe,Aye  and 
aperture  fields  EXe>Eye  of  each  identical  element  in  the  four-horn 
array.  Pattern  separability  and  ideal  element  excitation  are  addi¬ 
tional  necessary  assumptions.  Starting  with  the  general  expression 
for  the  radiation  field  of  the  array  in  wavenumber  coordinates. 


F  fit  v  \  =  jr  k  %ra  0j6(dk+dk)  jg  (-d  k  +d  k  ) 

-ffk  x’V  ^ffUx* VL  1  e  0  xx  yy  +  a2e  o  xx  yy 


(1) 


,  -JB  (d  k  +d  k  ),  -j8  (-d  k  +d  k  )n 

+  a,  e  o  x  x  yy+a,  eo  xx  yyj, 
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:  -1  -1-*  ‘ 


it  can  be  shown  that  the  H-plane  plane  wave  spectra  of  the  element 
may  be  found  from 


W0'  +  'Wk*-°> 


A  l\t  O'!  A  VL*r***  ^  _ _ 

xe'  x*  '  (1+k  )  (4  cosg  d  k  +2  sinf3  d  k  ) 


(2) 


O  X  X 


O  X  X 


„  .  .  W1”  +  'WV°> 

ye'  x*°  (l*fk  )(4  cos6  d  k  +2  sing  d  k  ) 


(3) 


O  X  X 


0  X  X 


where  the  numerators  are  the  measured  amplitude  data  and  where 
kx*sin6  in  the  H-plane.  Similar  expressions  hold  for  the  E-plane. 

The  above  expressions  for  A^e,Aye  are  valid  only  over  the  range 
^■±1;  i.e.,  0*»±9O°.  In  order  to  obtain  (via  Fourier  transform)  the 
corresponding  element  aperture  fields,  evanescent  energy  must  be 
added  to  the  plane  wave  spectra  according  to  well-define  criteria 
[3]  .  The  criteria  adopted  here  is  that  the  element  aperture  fields 
are  limited  in  spatial  extent  to  the  circular  region  occupied  by  the 
horn  aperture j  i.e.,  are  band-limited. 

The  first  estimate  to  the  field  Exnf  (x,o)  is  obtained  by  Fourier 
transforming  Axe(kx,0).  Energy  outside  |x|<a  is  truncated,  and  the 
resulting  field  is  inverse  transformed  to  yield  a  spectrum  A^Ck^O). 
The  energy  outside  j  k^ | <.1  is  taken  as  the  first  estimate  to  the  evanes¬ 
cent  energy  of  A^.  This  evanescent  energy  is  added  to  A^,  and  the 
process  is  repeated  until  a  very  high  percentage  of  energy  in  Exnf  is 
concentrated  in  the  region  |x|^a.  The  procedure  is  repeated  for  the 
other  field  component.  A  complete  element  aperture  field  can  be  ob¬ 
tained  to  within  a  multiplicative  constant  as  E  (k  ,k  )*E  (k  ,0) 

E  (0,k  ).  xe  x  y  xe  x 

ye  ’  y 

Results 


Figure  1  shows  a  comparison  of  the  true  and  synthesized  differ¬ 
ence  elevation  patterns  for  a  four-horn  array.  The  true  data  were 
generated  using  a  TEjjcircular  waveguide  mode  in  each  element.  The 
mode  was  rotated  5°from  the  vertical  to  produce  the  low  amplitude 
cross  (x)  component.  Concentration  factors  of  (1.-10'^)  and  (l.-10~^) 
were  obtained  for  the  primary  (y)  and  cross  (x)  components,  respec¬ 
tively,  to  produce  the  excellent  synthesis  results  shown. 
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FIGURE  1.  COMPARISON  OF  TRUE  AND  SYNTHESIZED  (')  E-PLANE  DIFFERENCE 
PATTERNS  OF  FOUR-HORN  MONOPULSE  ANTENNA. 
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It  is  well  known  that  the  receiving  and  transmitting  patterns  of 
an  antenna  in  free  space  are  identical.  In  what  follows,  it  is  theoreti¬ 
cally  demonstrated  that  the  receiving  and  transmitting  patterns  for  an 
antenna  enclosed  by  a  dielectric  radome  are  also  identical,  and  that  the 
receiving  and  transmitting  formulations  of  radome  analysis  [1]  lead  to 
identical  results. 

The  Huygens-Fresnel  principal  [2]  and  the  Lorentz  reciprocity 
theorem  [3]  are  used  to  establish  the  equivalence.  The  theory  is  then 
applied  to  show  that  the  effect  of  the  scattered  fields  on  receiving 
response  is  identically  zero.  A  method  is  suggested  whereby  discrepancies 
caused  by  an  unrealizable  edge  condition  for  an  aperture  illumination 
function  can  be  resolved  for  the  receiving  case  by  utilizing  results 
already  obtained  for  the  transmitting  case  [4] . 

Theory 

Consider  a  linear  antenna  operating  in  the  sinusoidal  steady  state 
(e‘,Wt  time  variations  suppressed)  and  enclosed  by  a  dielectric  radome. 

Let  a  closed  outer  surface  S  coincide  with  the  outer  surface  of  the 
radome  so  that  when  the  antenna  is  energized,  the  electric  field  at  any 
point  (x'jY'jZ')  outside  the  radome  is  given  by  [2]. 


E(x’,y’,z’)  =  —  J  [-jwyiHn  x  +  (n  x  E^)  x  Vf  (n’E^V^ds  (1) 
s 

where  are  the  fields  produced  on  S  by  the  antenna,  n  is  a  unit 

vector  normal  to  S  and  pointing  into  the  unbounded  region,  and  ip  is  the 
free-space  Green's  function 


<P  = 


(2) 


and  where  r  is  the  distance  from  the  (equivalent)  source  point  (x,y,z) 


on  S  to  the  point  of  observation  (x^y'.z'}  as  illustrated  in  Figure  1. 
For  the  directions  assigned  to  the  unit  vectors  r  and  R  in  Figure  1, 
the  expression  for  gradient  i/j  is  {2] 

l  e“jkr  ^ 

=  (jk  +  r  (3) 

where  k  =  2ir/X  and  where  X  is  the  wavelength.  The  constitutive  para¬ 
meters  of  the  homogeneous  unbounded  medium  are  z  and  p . 

After  Collin  [5] ,  let  an  infinitesimal  current  element  of  strength 
1^,  length  A£,  and  orientation  nfa  be  located  at  <x' ,y’ ,z')  .  Let  the 

fields  produced  on  S  by  the  current  element  be  given  by  E  ,H  .  An 

R 

application  of  the  reciprocity  theorem  to  the  unbounded  region  yields 


/(E  x  H  -  E  x  H  )  •  ndS  =  -  E  •  rrl  At 
— T  — R  — R  — T  —  b  b 


A  second  application  of  the  reciprocity  theorem  to  the  source-free 

volume  V'  inside  the  radome  enclosed  by  S  and  by  a  surface  S'  surrounding 

the  antenna  yields  the  result  that  the  complex  voltage  response  V  of 

R 

the  antenna  to  the  fields  E,H  ,  incident  on  S  is  proportional  to  the 

R  R 

integral  in  Equation  (4);  hence, 


vr  =  V^,nb  =  c/<£ r  *  Hr  -  ER  X  .nds 


where  E_  is  given  by  Equation  (1)  in  the  absence  of  the  current  element 

and  where  C  is  a  complex  constant. 

When  the  point  of  observation  is  a  great  distance  R  from  the  origin 

0  of  source  coordinates,  the  fields  E , H  produced  on  S  by  the  current 

R  R  ^ 

element  are  those  of  a  plane  wave  propagating  in  the  -R  =  k  direction, 

a 

and  is  given  by 

E  x  R  k  x  E 
„  _  ~ R  _  a  — R 


where  g  is  the  intrinsic  impedance  of  the  unbounded  medium.  The  constant 

C  in  Equation  (5)  can  be  chosen  so  that  E  is  given  by 

R 


19 


(7) 


e-jkR  _  jkka.r' 

— r  =  ^  —T~  V 


where  £'  is  the  position  vector  from  the  origin  0  to  the  source 
coordinates  (x,y,z)  as  shown  in  Figure  1.  The  magnetic  field  follows 
from  Equations  (6)  and  (7)  as 


-jkR  jkk  £' 

Hr  =  jk  C  <nb  x  R) 


(8) 


where  the  relation  k  =  wp/n  has  been  used. 

The  next  step  in  the  development  consists  of  evaluating  Equation 
(5)  using  the  expression  for  E  from  Equation  (1) ,  under  the  condition 
of  large  R;  i.e.,  far-field  approximation.  When  this  is  done,  the 
response  of  the  antenna  to  the  incident  plane  wave  is  given  by 

lAl 
b 

R  4tt  •£  1  R  “  "b 


V  = 


b —  /*  e~jkR  3kka*r' „ 

J  [-ju>p  — - —  e  n^  •  ( n  x  +  n,_»(n  x  H£,) 


-T 


"  iicP_jkR  jkk  •£' 

x  R  - —  e  +  n  (n-E  )  ]  dS 

R  b  ff  — T 


(9) 


Where  the  following  asymptotic  relations  have  been  used: 

r~R  -  R.r' 


(10) 


-jkr  -  jkR  jkk  -r' 

e  J  e  „  a  — 

-  ~  — - —  e 


(ID 


-jkR  jkk  •£' 
Vi =  Rjk  r 


(12) 


Using  the  vector  identify  £*bxc  =  b*£X£,  the  first  term 
T^  in  the  integrand  of  Equation  (9)  can  be  written  as 


-jkR 


T1  "  — 


n  •  H_  x  n^e 


jkk  .£' 

3. 


(13) 


Comparison  to  Equation  (7)  reveals  that 


T  =  E  x  H  •  n 
1  -R  -T 


(14) 
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i.e.,  T  is  identical  to  the  second  term  in  the  integrand  of  Equation 
(5)  . 

Repeated  application  of  the  same  vector  identity  to  the  second 
term  in  the  integrand  of  Equation  (9)  yields 

t2  =  n  *  [*Lr  x  Wff  x  n^)]  (15) 

or 

A 

^  *■  j3cir  jkk  •  2T 1  ^ 

T2  =  n  ■  [^r  x  (jk  e  a_  Rxr^)]  (16) 


Comparison  of  the  expression  in  parenthesis  in  Equation  (16)  with  Equation 
(8)  shows  that  is  equal  to  the  negative  of  the  first  term  in  Equation 
(5) .  Collecting  results  yields  the  desired  final  result  as 

-  A  A 

V  =  c/(E  xH  -  E  xH).  ndS  =  n  •  E 
R  J  —T  — R  — R  — T  b  — 

S 

A  A  A  A 

=  n bJ  [- jwpijj (n  x  i^)  +  (n  x  E^,)  x  Vipff  +  (n-E)ViHds  (17) 

S 

for  the  case  of  large  R  (far-field  approximation) .  It  is  noted  that 
the  third  term  in  Equation  (9)  is  identically  zero  since  n^  is  perpen¬ 
dicular  to 

The  result  in  Equation  (17)  establishes  that  the  receiving  and 
transmitting  patterns  of  the  antenna/radome  combination  are  identical. 
Furthermore,  it  establishes  the  equivalence  of  the  receiving  and  trans¬ 
mitting  formulations  of  radome  analysis  under  far-field  conditions.  The 
analysis  also  shows  the  conditions  of  equivalence;  viz.,  that  the 
fields  on  the  outer  surface  of  the  radome  must  be  known.  In  this  respect, 
the  receiving  formulation  may  offer  some  advantages  over  the  transmitting 
formulation  since  it  does  not  depend  on  the  existence  of  a  homogeneous 
unbounded  medium  outside  the  surface  S . 

Application 

A  concise  application  of  the  foregoing  analysis  is  to  show  that 
the  response  of  the  antenna  to  a  plane  wave  depends  on  only  the 
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Examination 


i  1  S  S 

incident  field  E,.  H  and  not  on  the  scattered  fields  E  ,  H  . 

— Ft  -R  — R  — R 

of  Equation  (9)  shows  directly  that  only  the  incident  plane  wave  is 
involved.  This  is,  of  course,  in  keeping  with  the  reciprocity  theorem 
in  that 


/(E  x  H?  -  E®  x  H  ) • ndS  HO  (18) 

— T  — R  — R  — T 

S 

when  the  source  free  volume  V  is  considered. 

A  second  application  of  the  analysis  concerns  the  problem  of 
resolving  differences  in  responses  for  aperture  distributions  having 
bounded  support  by  using  line  charge  distributions.  For  the  trans¬ 
mitting  case,  the  method  of  Stratton  and  Chu  [4]  can  be  readily  utilized. 
Moreover,  the  same  method  can  be  applied  in  the  receiving  case  by  utilizing 
the  above  analysis  results  as  will  be  done  in  a  later  paper. 
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Chapter  1 


INTRODUCTION  AND  SUMMARY 


1-1.  Introduction 

This  report  documents  a  ray  tracing  radome  analysis  computer  program 
written  in  Fortran  IV  for  use  on  the  IBM  3033  computing  system  at  Johns 
Hopkins  University  Applied  Physics  Laboratory.  The  program  was  developed 
at  Georgia  Institute  of  Technology  over  the  past  four  years;  however, 
considerable  development  work  in  computer  aided  radome  analysis  has  taken 
place  here  prior  to  that  time  [1-7] . 

This  analysis  package  is  intended  to  support  the  on-going  radome 
technology  program  at  APL  and  is  documented  herein  for  that  purpose; 
however,  it  is  important  to  point  out  that  this  effort  was  done  in  con¬ 
junction  with  the  on-going  radome  analysis  work  at  Georgia  Tech  under 
grant  AFOSR-77-3469 [8] .  It  is  anticipated  that  this  report  will  be 
available  to  other  government  contractors  through  the  National  Technical 
Information  Service  as  part  of  the  technology  transfer  functions  of  the 
supporting  agencies. 

The  report  is  organized  by  chapters  according  to  the  approximate 
order  in  which  the  subprograms  are  called,  and  each  chapter  describes 
one  subprogram.  Each  chaper  is  essentially  self-contained  since  it  is 
meant  to  serve  as  the  complete  documentation  on  a  single  subroutine. 
References  are  provided  at  the  end  of  each  chapter.  In  some  cases, 
figures  are  duplicated  in  different  chapters  for  completeness.  Each 
chapter  is  terminated  with  the  listing  of  the  subroutine. 


Chapter  2  describes  the  main  program  and  instructions  for  its  use. 
Chapters  3  through  28  describe  the  thirty  four  subroutines  required  for 
execution,  including  those  for  producing  Calcomp  pattern  plots  and  three- 
dimensional  plots.  Appendices  A  through  D  present  computed  results  for 
four  test  cases  for  use  in  verifying  correct  operation  on  other  systems. 
These  results  were  obtained  on  the  Cyber  70/74  computing  system  at  Georgia 
Tech.  The  remaining  part  of  this  chapter  describes  background  of  the 
program  development  and  summarizes  the  features  of  the  computer  analysis . 

Two  companion  reports  are  also  of  interest  [9,  10] .  In  the  first, 
a  method  of  analysis  based  on  the  Huygens-Fresnel  principle  (surface 
integration)  is  described  which  makes  use  of  most  of  the  computer  soft¬ 
ware  described  in  the  present  report.  The  second  report  compares  the 
results  obtained  using  these  two  computer-aided  analyses  with  those 
obtained  experimentally  at  35  GHz  for  various  combinations  of  sizes  of 
monopulse  antennas  and  tangent  ogive  radomes  [8] .  The  latter  report 
represents  the  culmination  of  this  radome  analysis  research, 

1-2 .  Background 

Development  of  the  radome  analysis  computer  program  (RACP)  was 
i  ii  i  L  l.<  Led  in  in  an  effort  to  include  the  dice Cs  of  the  rudoinc  on 

a  ground  mapping  radar  [1] .  A  three-dimension  geometry  and  vector  field 
formulation  were  used.  A  plane  wave  Si>octruin  (l’WS)  representation  of 
the  radiation  from  the  antenna  greatly  facilitated  the  computations 
since  the  Fast  Fourier  Transform  (FFT)  could  be  used.  The  program  was 
used  to  compute  power  patterns  on  the  ground  for  many  different  cases 
of  antenna/missile  orientations.  From  these  data,  the  effects  of  the 
radome  on  pattern  shape ,  power  loss  and  VSWR  were  determined. 
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Monopulse  tracking  antennas  were  next  introduced  into  the  computer 
analysis  to  evaluate  radome  materials  and  shapes  for  seeker  systems  in 
the  8-18  GHz  band  [2] .  Tangent  ogive  shapes  of  various  fineness  ratios 
were  analyzed.  Monolithic  and  multilayer  wall  structures  were  used. 
Algorithms  were  developed  to  compute  boresight  errors  from  the  sampled 
data  difference  patterns  in  two  orthogonal  planes.  A  modification  of  this 
program  was  also  used  to  conduct  a  trade-off  and  development  study  for 
the  Multipurpose  Missile  (MPM) ,  later  known  as  ASALM  [3] . 

The  next  step  in  the  development  of  RACP  came  in  1977  with  the 
introduction  of  a  conical  scan  tracking  antenna  into  the  analysis  [4] . 

This  antenna  necessitated  a  reformulation  of  the  analysis  from  the  trans¬ 
mitting  formulation  used  earlier  to  a  receiving  formulation.  The  big 
advantage  offered  by  the  latter  is  that  the  antenna  response  can  be 
calculated  for  only  one  direction  of  arrival  of  the  target  return  (plane 
wave)  .  In  the  former,  the  FFT  automatically  computes  ’’responses"  for 
many  directions  of  arrival  and,  hence,  is  computationally  slower. 
Subsequent  versions  of  the  program  have  used  the  same  receiving  for¬ 
mulation  with  monopulse  and  other  types  of  antenna  models. 

The  computed  results  obtained  with  the  receiving  and  LransmiLting 
formulations  are  not  always  the  same  15] .  A  computed-aided  analysis 
which  utilizes  the  Huygens-Fresncl  principle  [6,  7]  is  generally  con¬ 
sidered  to  be  more  accurate  than  the  two  methods  already  mentioned,  but 
requires  considerably  more  computation  time  that  may  not  be  warranted 
in  all  cases.  A  research  program  is  now  underway  at  Georgia  Tech  whose 
objective  is  to  establish  the  accuracies  and  ranges  of  validity  of  those 
three  methods  of  radome  analysis  [5]. 
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1-3. 


Description  of  the  Analysis 


The  current  version  of  the  ray  tracing  analysis  computer  program 
utilizes  a  receiving  formulation  based  on  the  Lorentz  reciprocity  theorem 
15].  A  plane  wave  of  selectable  linear  or  circular  polarization  is 
assumed  incident  on  the  outside  of  the  radome  and  is  represented  by  a 
system  of  parallel  rays.  There  is  one  ray  for  each  sample  data  point  in 
the  antenna  aperture  inside  the  radome.  Each  ray  is  traced  from  the  point 
where  it  impinges  on  the  outside  surface  to  the  corresponding  aperture 
point.  The  electric  and  magnetic  fields  E^,  fr  associated  with  each  ray 
are  weighted  by  the  flat  panel  transmission  coefficients  T^,  T ^  as  de¬ 
termined  by  the  unit  normal  n,  the  direction  of  propagation  k,  and  the 
dielectric  properties  of  the  radome  wall.  The  weighted  incident  fields 
E|,  at  each  aperture  point  are  then  used  in  the  following  integral  to 
obtain  the  complex  voltage  response  V^_  of  the  antenna  as 
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where  E^,  are  the  aperture  fields  when  the  antenna  is  transmitting, 

C  is  a  complex  constant,  and  z  is  the  unit  vector  normal  to  the  xy 
(aperture)  plane.  For  digital  computer  implementation,  the  integral  in 
Equation  (1)  reduces  to  a  double  summation,  and  the  equal-area  elements 
dxdy  Ijccome  AxAy  and  can  be  absorbed  into  the  constant  C. 

In  its  present  form,  the  program  accommodates  only  one  radome 
shape;  viz.,  the  tangent  ogive.  The  length,  diameter  and  fineness  ratio 
are,  of  course,  all  variable  in  the  input  data.  Monolithic  and  multi¬ 
layer  wall  configurations  can  be  analyzed;  however,  only  uniform  wall 
configurations  whose  properties  do  not  vary  from  point  to  i>oint  on  the 


wall  can  be  handled.  Provisions  are  made  to  allow  for  a  metal  tip  on 
the  radome  whose  effect  is  aperture  blockage. 

The  geometry  subroutines  provide  for  three  separate  coordinate 
systems  and  the  point  and  vector  transformations  among  them.  A  reference 
coordinate  system  is  provided  to  orient  the  antenna/radome  combination 
with  respect  to  other  bodies.  The  coordinate  systems  for  the  antenna 
and  the  radome  comprise  the  other  two  systems.  Boresight  error  and 
pattern  computations  are  carried  out  and  exjjrersed  in  the  antenna  coor¬ 
dinate  system. 

The  primary  outputs  of  the  program  are  boresight  error  (mrad.), 
boresight  error  slope  (deg. /deg.),  gain  loss,  and  when  selected,  prin¬ 
cipal  plane  patterns.  Outputs  include  both  printing  and  plotting 
(Calcomp) .  Plotting  options  allow  for  selection  of  aperture  fields 
with  and  without  the  radome.  A  feature  is  also  provided  to  either  ob¬ 
tain  or  suppress  intermediate  calculated  results  for  debugging  purposes. 

Boresight  error  calculations  for  monopulse  antennas  are  carried 
out  by  setting  the  first  target  return  at  a  known  direction  within  a 
few  degrees  of  true  boresight.  The  responses  in  the  two  difference 
channels  and  the  sum  channel  are  then  computed  and  stored.  Another  set 
of  responses  for  a  return  180°  away  from  the  first  is  computed  next. 

The  two  sets  of  data  are  then  used  to  construct  a  linear  tracking  model 
in  the  two  orthogonal  planes,  and  the  process  is  repeated  until  a  bore¬ 
sight  null  is  indicated.  The  true  direction  of  arrival  of  the  plane 
wave  at  this  point  represents  the  boresight  error  directly. 

The  current  subroutine  used  to  characterize  the  antenna  permits 
selection  of  various  polarizations  and  two  aperture  distributions.  A 
uniform,  circular  aperture  distribution  having  vertical,  horizontal  or 


circular  (LHR  or  RHC)  polarizations  is  one  combination.  The  second 
distribution  is  a  tapered  rectangular  distribution  having  vertical  polar¬ 
ization  as  found  in  flat  plate  antennas .  This  basic  subroutine  would  not 
be  difficult  to  modify  to  accommodate  other  distributions,  such  as  rec¬ 
tangular  aperture  with  cosine  taper. 

Computation  time  is  independent  of  radome  size  but  depends  on  the 
number  of  samples  used  in  the  aperture.  For  256  sample  points  (16  X  16 
array) ,  the  time  to  compute  the  received  voltages  in  the  three  channels 
is  1.5  seconds. 

The  program  is  organized  as  a  main  program  and  a  number  of  sup¬ 
porting  subroutines,  all  written  in  Fortran  IV.  The  complete  program, 
including  plotting  software,  contains  thirty  four  subroutines.  The  core 
storage  required  for  the  complete  program,  including  all  library  and 
system  I/O  routines,  is  just  over  46,000  (decimal)  words.  Integer,  real 
and  complex  variables  and  arrays  are  utilized.  Single,  double  and  three- 
dimensional  data  arrays  are  present.  Only  single  precision  variables  and 
computations  are  required  with  the  60-bit  word  available  on  the  Cyber  70 
at  Georgia  Tech. 
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3-0  3.  Prerequisites:  Graduate  standing.  Maxwell's  equations,  antennas, 

vector  analysis,  and  multivariable  calculus. 

Radome  analysis  is  the  application  of  basic  electromagnetic  theory  to 
determine  the  effects  of  protective  dielectric  structures  on  the  electrical 
characteristics  of  antennas  so  enclosed. 

Course  Objective:  To  enable  the  participants  to  derive  selected  important 
relations  from  Maxwell 1 s  equations  and  to  apply  them  to  the  solution  of 
practical  electromagnetic  problems  dealing  with  the  interaction  of  radiating 
structures  and  their  dielectric  enclosures;  viz.,  radomes  and  electromagnetic 
windows . 

Text:  (1)  S.  Silver  (ed.).  Microwave  Antenna  Theory  and  Design,  Volume  12 

of  MIT  Radiation  Laboratory  Series  (printed  by  Georgia  Tech 
Bookstore) 

(2)  Selected  journal  papers  on  reserve  at  library 

Topical  Outline 

Definition  of  the  Problem 

History  of  radome  development 
Electromagnetic  windows 
Antenna  systems 

Aerodynamic  and  mechanical  considerations 

Survey  of  Existing  Analysis  Methods 

Geometrical  optics  (Silver) 

Equivalent  aperture  (Kilkoyne) 

Integral  equation  (Van  Doeren) 

Equivalence  and  -'nduction  approach  (Tricoles) 

Corrections  for  wall  curvature  (Tavis) 

Huygens  principle  (Paris) 

Plane  Wave  Spectrum  (Wu  and  Rudduck) 

PWS  and  Equivalent  Aperture  (Joy  and  Huddleston) 

Spherical  Wave  Expansions  (Chesnut) 

Lorentz  reciprocity  (Huddleston  and  Joy) 

General  Theory  of  Radome  Analysis 

Transmitting  formulation 
Receiving  formulation 
Integral  equation 
Equivalence  of  formulations 
Electromagnetic  windows 
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Electromagnetic  Theory 

Vector  Green’s  theorem 
Kirchof f-Huygens  principle 
Reciprocity  theorem 
Equivalence  theorems 
Geometrical  optics 
Scattering 

Plane  waves  and  infinite  dielectric  sheets 
Modal  expansions  for  antenna  radiation 

Special  Considerations 

Coordinate  systems  and  transformations 
Intersections  of  rays  with  surfaces  of  revolution 
Surface  integration 
Approximations 
Comparisons  of  results 

Electromagnetic  Windows 

Dielectric  slab-covered  waveguide 
Integral  equation  formulation 
Method  of  moments  solution 
Radiation  and  impedance 

Boundary  Value  Problems 

Dipole  within  concentric  spherical  shell 
Spherical  antenna  with  dielectric  sheath 
Radiation  from  conical  structures 

Ellipsoidal  antenna  with  confocal  dielectric  sheath 


